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ABSTRACT 

We present a two-dimensional grid-based hydrodynamic simulation of a thin, viscous, locally-isothermal 
corotating disk orbiting an equal-mass Newtonian binary point mass on a fixed circular orbit. We study the 
structure of the disk after multiple viscous times. The binary maintains a central hole in the viscously-relaxed 
disk with radius equal to about twice the binary semimajor axis. Disk surface density within the hole is reduced 
by orders of magnitude relative to the density in the disk bulk. The inner truncation of the disk resembles the 
clearing of a gap in a protoplanetary disk. An initially circular disk becomes elliptical and then eccentric. Dis- 
turbances in the disk contain a component that is stationary in the rotating frame in which the binary is at rest; 
this component is a two-armed spiral density wave. We measure the distribution of the binary torque in the disk 
and find that the strongest positive torque is exerted inside the central low-density hole. We make connection 
with the linear theory of disk forcing at outer Lindblad resonances (OLRs) and find that the measured torque 
density distribution is consistent with forcing at the 3:2 (m = 2) OLR, well within the central hole. We also 
measure the time dependence of the rate at which gas accretes across the hole and find quasi-periodic structure. 
We discuss implications for variability and detection of active galactic nuclei containing a binary massive black 
hole. 

Subject headings: accretion, accretion disks — binaries: general — black hole physics — galaxies: nuclei — 
hydrodynamics 



1. INTRODUCTION 

Circumbinary disks have been observed in young stellar 
binaries, such as the spa tially-resolved pre-main-sequence 
binary star GGTau (e.g., iMcCabe. Duchene. & Ghezl 120021: 
iKrist et aTl 120051) . Circumbinary disks should also form 
around binary massive black holes when interstellar gas ac- 
cretes into the binary's dynamical sphere of influence (e.g., 
Artvmowicz & Lubow 1996; Ivanov, Paoaloizou, & Polnarev 
1999 1 ; lArmitage & Nataraianll2002t IMilosavlievic & Phinnevl 
2005b . It is generally expected that the dynamics of a binary 
embedded in an accretion disk is similar to the dynamics of 
a planet embedded in a protoplanetary disk. The planet or- 
bits the central star while exchanging gravitational torques 
with, and accreting from, a gaseous circumstellar disk. If 
the radius of the planet's Roche lobe exceeds the pressure 
scale height of the disk, the planet clears a gap in the disk 
(e.g.jGoldreich & Tremainel[l9 80; Takeuchi, Mivama. & Linl 
ll996tlLubow. Seibert. & Artymowiczl 19991) . When a similar- 
mass binary is surrounded by a thin corotating disk, this gap- 
opening condition is automatically satisfied. 

An understanding of circumbinary accretion flows is es- 
sential for future identification of binary massive black holes 
in astronomical surveys. Massive binary black holes have 
only tentatively been identified in astronomic al observations 
( Komossa 2003, 120061 : iRodriguez et al.l l2006 report evidence 
for a binary with projected separation of 7.3 pc), but the 
theoretical case for their existence is very strong. A mas- 
sive binary black hole forms during the conclusion of a 
galaxy merger and persists at parsec or subparsec separations 
over a cosmological time. Under propitious circumstances, 
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gravitational radiation ultimately induces coalescence be- 
tween the black hol es (Begelman. Blandford. & Reesl 119801 
Merritt & Milosavljevic 2005, and references therein). 

At parsec separations, the black holes can be optically re- 
solved only in nearby galaxies. Better resolution can be 
achieved at radio frequencies, but binary radio cores can be 
confused with structures in radio jets and chance superposi- 
tions. An unresolved binary can be detected by associated 
characteristic spectral and time-dependent signatures if a cir- 
cumbinary accretion flow is present. For example, asymme- 
tries in the accretion flow and time-dependent accretion onto 
the component black holes are sources of variability that can 
be detected in monitoring surveys such as with the Large Syn- 
optic Survey Telescope (LSST) 6 . 

To predict the spectral signatures of massive binary black 
holes, one must analyze the hydrodynamic structure of the 
circumbinary accretion flow and then synthesize the radia- 
tive spectra emitted by the flow. In this work we carry out 
the hydrodynamic analysis in the case of a thin corotating 
disk . We defer the spectral synthesis to subsequent work (se e 
also IMilosavlievic & Phinnevll2~005l; iBogdanovic et al.ll2007l) . 
The search for massive binary black hole spectral transients 
is most effective if one can identify robust transients that re- 
peat continuously over an extended period. Transients that are 
restricted to a single moment in the life of a binary, such as 
its formation, will be rare and difficult to interpret. Quasi- 
steady-state circumbinary accretion disks are particularly at- 
tractive as they may persist long enough to render their astro- 
nomical detection probable. 

Previous numerical simulations of circumbinary disks 
based on the smooth e d part icl e hydrodynamics m ethod 
dArtvmowicz & Lubowl 119941; lEscala et al.l 120051) and 
grid-based hydrodyna mics with explicit bulk viscosity 
( Gunther & Kleyl [20021) suggest that the surface density of 
a circumbinary disk drops sharply inward of r ~ 1.5- 2a, 

6 http://www.lsst.org 
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where a is the binary semimajor axis. The binary resides 
within a low-density hole, also called the circumbinary gap, 
at the center of the disk. Although the gas density within 
the hole is low, some mass transfer does take place from th e 
disk onto the point masses (iGunther. Schafer. & Klevll200l . 
However the accretion rate is below the rate that would be 
expected neglecting the torque that the binary exerts on the 
disk. The reduced accretion was detected in simulations 
of massive planets embedded i n prot o planetary disks (e.g, 
Lubow. Seibert. & Artvmowiczl Il999t iLubow & D'Angelol 
2006). 

The simulations of disks with planets also indicate that an 
initially circular disk becomes eccentric over time, even if the 
orbit of the planet is circula r (Papaloiz ou. Nelson. & M asset 
l200UlKlev & Dirksenll2006l) . The growth of di sk eccentricity 
has been attributed to a dynam ical instability (Lubow| ["l991t 
ID'Angelo. Lubow. & Bate(l2006l) driven by the planet's tidal 
potential. In turn, the disk eccentricity can excite eccen- 
tricity in the binary, since the two ec centricities are coupled 
(Papaloizo u. Nelson. & M asset 2001). Even a weak resid- 
ual eccentricity can be detected in the gravitational radiation 
emitt ed during black hole coalescence ( Armitage & Nataraianl 
2005; see also Arty mowicz et al.ll 19911) . Therefore an under- 
standing of the circumbinary disk dynamics is key for elec- 
tromagnetic and gravitational wave observations of coalesc- 
ing massive black holes with the Laser Interferometer Space 
Antenna {LISA). 7 Simultaneous electromagnetic and gravi- 
tational detection can be used to measure the geometry of 
the universe and constrain properties o f dark energy (e.g., 
iHolz & Hughesll2005tlKocsis et al.ll2006T) . 

We present high resolution simulations of a corotating vis- 
cous circumbinary disk around a circular Newtonian equal- 
mass binary. The disk lies in the orbital plane of the binary. 
This is an astrophysically-motivated simplifying assumption; 
alignment of the orbital plane of the inner disk and that of 
the binary will take place because differential precession of a 
non-aligned disk caused by the binary's mass quadrupole mo- 
ment forces the inner par t of the disk to align with the orbital 
plane of the binary (e.g jLarwood & Papaloizou! 19971) . 8 

In § |2] we describe the computational method and the 
binary-disk setup. In § [5] we present our results. We study 
the disk eccentricity, the spatial distribution of torque depo- 
sition in the disk, and the mass accretion across the central 
hole in the disk. In § |4] we compare the torque distribution 
measured in the simulation with predictions of the theory of 
linear response of tidally forced disks. We also present a ten- 
tative interpretation of the eccentricity excitation in terms of a 
mean motion resonance. In §|5]we discuss implications for as- 
tronomical detection of binary massive black hole mergers in 
progress. Finally, in §|6]we summarize the main conclusions. 



2. THE CODE AND INITIAL CONDITIONS 

The simulation w as carried out with the publicly available 
code ASC FLASH dFrvxell et al.ll2000l) in two spatial dimen- 
sions. The mesh refinement capability of FLASH allowed us 
to use high grid resolution in the central region of the disk 
(within a radius equal to 10 binary semimajor axes) and lower 
resolution in the outer disk. We implemented a static grid 
with constant resolution in the cf) direction. The piecewise- 

7 http://lisa.nasa.gov 

8 The orbital plane of the binary will itself precess due to the relativistic 
spin-orbit coupling, if the binary components have spins misaligned with the 
binary's angular momentum. 



parabolic method module implemented in FLASH is ideally 
suited for treating supersonic dynamics in the vicinity of the 
binary. The simulation was carried out in polar (r,(j>) coordi- 
nates in the inertial frame. 

We implemented several small changes in FLASH. We in- 
troduced a Newtonian binary point particle of total mass M 
on a fixed circular orbit centered at the origin. The gravity of 
the binary acting on the disk fluid is evaluated directly; the 
gravity of the fluid is not calculated. The integration domain 
is a circular annulus with an inner boundary at r m ; n = a and 
outer boundary at r max = 100a. The binary point masses orbit 
at the distance a /2 from the center; therefore the point masses 
move entirely within the excised circular region r < r m ; n . An 
obvious limitation implied by the excision is that we do not 
study the accretion of gas onto the individual point masses, 
but we do measure the net rate of mass transfer across the 
inner boundary. 

The fluid was allowed to flow freely out of the integration 
domain at the inner and outer boundary, but the flow into the 
domain was forbidden. Specifically, the boundary conditions 
at r m ; n and r max were such that the radial derivative of the sur- 
face density, <9X/<9r, and of the azimuthal velocity, dv^/dr, 
were set to vanish at each boundary. The boundary condition 
imposed on the radial velocity v r was a function of the direc- 
tion of the flow at the boundary. For outflow, v,(r m ; n ) < 
or v r (>"max) > 0, the radial derivative of the radial velocity, 
dv r / dr, was set to vanish; for inflow, the value of v r in the 
boundary ("ghost") zones just adjacent to r m ; n and r max was 
set to vanish, thus preventing inflow. 

We implemented a locally isothermal equation of state in 
which the sound speed is a fixed function of radius c s (r). 
The disk has a constant pressure scale height to radius ratio 
h/r = Cs/rQtep = 0.1, where Sl kep = (GM/r 3 ) 1 / 2 is the Kep- 
lerian frequency in the disk; therefore the sound speed de- 
creases with radius c s oc r -1 ' 2 . Since this equation of state is 
applied exactly, the total energy of the fluid is not conserved 
and all shocks can be considered radiatively isothermal. 

We also implemented an explicit viscous shear stress term 
in the equation governing the transport of the 0-momentum. 
In doing so, we kept only the (r, (j>) component of the stress 
tensor, since this is the component that communicates viscous 
torque in the disk and is responsible for accretion. The corre- 
sponding vertically-integrated Navier-Stokes equation reads 

dt r dr r 86 rd<b 



1 



dP | d(r 2 a r ^ 
rd<j) r 2 dr 



(1) 



where £ is the surface density in the disk, 



$(r,0) = 



\GM 



[r 2 +\a 2 -ra cos(6> - W)] 1/2 



(2) 

[r 2 +ifl 2 + racos(6'-r2 bin f)] 1 / 2 

is the Newtonian gravitational potential of the binary point 
mass, fibin = (GM/fl 3 ) 1 / 2 is its angular frequency, P = Ec 2 
is the vertically-integrated pressure, and v r and are the ra- 
dial and the tangential component of the fluid velocity, respec- 
tively. 

The (r, 6) component of the stress tensor equals 

4(-) +i Sl< ® 

or \ r / r do 
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FIG. 1. — Surface density of the disk a function of radius in units of the 
binary semimajor axis. The curve peaking on the right is the initial disk 
density profile; the one peaking on the left is density profile after ~ 4,000 
binary orbits averaged over ~ 50 orbits. 



where v is the kinematic viscosity. On the right-hand side of 
equation ([3), the second term involving dv r /d(f> is normally 
much smaller than the first term involving diy^ /f)J dip. We 
drop the second term in the simulation. 

The viscosity is calculated based on the a-disk prescription 
v(r) = ac s (r) 2 /fikep(r) with a = 0.01. The use of the Keple- 
rian frequency in the evaluation of the viscosity is simply a 
prescription; the binary's mass quadrupole moment and posi- 
tive radial pressure gradients make the true angular frequency 
near the inner edge of the disk super- Keplerian 



kep 



3 a 2 
1+ 16 72 



1 dp 

r£ dr 



(4) 



Numerical resolution in the simulation was chosen to en- 
sure an accurate treatment of pressure waves excited in the 
disk. The binary's tidal potential can excite density waves 
at circular and eccentric outer Lindblad resonances (OLRs). 
These waves have wavelengths equal to a rational multiple of 
~ 27rc s (r)/i7bin ~ 0.6 a(r/a)~ l l 2 . The radial resolution was 
chosen such that the wavelength of a density wave is resolved 
by many resolution elements in the inner disk (within a few 
binary separations from the inner edge). Specifically, the 
simulation was carried on time-independent (non-adaptive) 
nonuniform meshes utilizing the PARAMESH adaptive mesh 
refinement package in FLASH. The spatial resolution was 
(Ar w 0.039 a, A0 w 0.0078 x 2?r) near the inner bound- 
ary (r ~ r m ; n ) and eight times lower (Ar « 0.31 a, A<j> » 
0.062 x 2tt) near the outer boundary (r ~ r max ). The resolution 
was uniform in the azimuthal direction at any fixed radius. 

The initial disk has radial surface density profile 



£(r)=£ 



exp 



(5) 



where r s = 10a, 8 = 3, £ = 2, and Eo is an arbitrary constant. 
The initial density profile is shown in Figure Q] The surface 
density decreases to a small value near the boundaries of the 
integration domain at r m i n = a and r max = 100a. 

The initial angular velocity of the fluid is corrected for 
the radial pressure gradient (see equation HJ and the ini- 
tial radial velocity is corrected for the viscous drift v r w 
2(r 2 VL'E)~ x d(r 2 ij r( f ) )/ dr. The initial velocity of the disk cor- 
rectly represents the quasi-steady state of an axisymmetric 



disk to the first order in h/r and a, and no significant tran- 
sients are observed at the beginning of the simulation. 

Since the inviscid hydrodynamical compone nt of FLASH 
has been extensively tested by other authors (iFrvxell et al.l 
2000), we are here mainly concerned with the correctness and 
stability of our implementation of explicit viscosity. To test 
the implementation, we simulated an identical initial accre- 
tion disk surrounding a single point particle of mass M at the 
origin. This disk was simulated at the same resolution as the 
circumbinary disk. The test disk remains axisymmetric and 
its evolution follows that prescribed by equation (fTJ on a vis- 
cous timescale, thereby validating our implementation of the 
explicit viscosity. 

Maintaining exact axisymmetry for many thousands of or- 
bits, as required for simulations lasting many viscous times, 
required care with the numerical grid. We set the FLASH 
parameter unbiased_geometry to clean the last bits of 
floating point numbers defining the numerical mesh. This was 
necessary because by default, due to numerical roundoff and 
truncation errors, the mesh spacing in the azimuthal direc- 
tion is not precisely uniform. This nonuniformity can seed 
nonaxisymmetric perturbations on many orbital timescales; 
unbiased_geometry enforces exact uniformity in the az- 
imuthal mesh spacing, thus eliminating spurious loss of ax- 
isymmetry which developed in test simulations run for a sin- 
gle point mass with and without viscosity. 

3. RESULTS 

3.1. Surface Density Profile 

The initial disk given in equation (0 is smoothly truncated 
at r ~ 3a. The binary tide experienced by the initial disk is 
weak and the disk initially accretes inward on a viscous time 
scale 

-i -l 



3 v 



- n 



, r \ 3/2 

1.MHH-) fbin, 



(6) 



where fbin = 27r/f2 bin is the orbital period of the binary. Near 
the inner edge the disk responds significantly faster to the vis- 
cous torque because of the presence of a steep density gra- 
dient. As the disk accretes inward toward the binary, the 
binary's tidal forcing of the inner edge of the disk strength- 
ens. When the viscously spreading disk arrives at a radius 
r ~ 2a, strong nonlinear disturbances develop in the circumbi- 
nary flow; we discuss these disturbances in detail below. 

The disk then settles into a quasi-steady state. Figure [T] 
shows the initial surface density profile and the time-averaged 
surface density profile after 4000 orbits. 9 The evolved disk 
is sharply truncated inward of r ~ 2a and the central region 
r <2a contains very little fluid. The surface density in this 
region is an exponential function of radius, 



(S)(r < 2a) w E exp(-15.2 + 5.95 r/a), 



where 



1 

27 



(7) 



(8) 



denotes the azimuthal average. While the azimuthally aver- 
aged inner density is low, accretion still occurs in streams 
of gas (Fig. |2|i which have radially infalling velocity compo- 
nents. 

9 Here and elsewhere in the text time averaging is carried out over ~ 50 
binary orbits. 
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FIG. 2. — Disk surface density at 1000 binary revolutions (top panel) and 
at 4000 revolutions (bottom panel); the black dots indicate the location of the 
point masses. The plotted region spans xja — ±5 and y/a = ±5, showing the 
inner 20% (in radius) of the simulation region. 



3.2. Eccentricity and Ellipticity of the Disk 

Dynamics of the inner disk is strongly nonlinear and time- 
dependent. The disk first becomes elliptical (Fig. [2] top) and 
then eccentric (Fig. [2] bottom). The elliptical state is sus- 
tained only on a single viscous timescale, but the subsequent 
eccentric state is persistent and apparently represents a quasi- 
steady state of the circumbinary disk. Therefore, we focus 
on the eccentric disk in the following discussion. We mea- 
sure the magnitude of the eccentricity, denoted by e\, and the 
disk ellipticity, denoted by e%, by taking a time average of the 
density-weighted m = 1 and m = 2 components, respectively, 
of the radial velocity, normalized by the density-weighted az- 
imuthal velocity 
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FIG. 3. — Time-averaged disk eccentricity ei (solid line) and ellipticity £2 
(dotted line) as a function of radius after 4,000 binary orbits. The two added 
in quadrature (e\ -t-el 1 )^ 2 are also plotted (dot-dashed line). Dashed thin line 
is the function exp(— r/a). 

nonlinear in the sense that the associated epicyclic motion is 
(mildly) supersonic. 

The disk eccentricity precesses slowly in the prograde di- 
rection, and the precession appears to be rigid at radii 2a < 
r <4a . The average precession rate of the eccentricity vector 
(the line of apsides) at these radii is 10 



fi prec = j arg (£v r e' ) 



;3.05 x 10" J O 



bin • 



(10) 



This can be compared with the apsidal precession rate 03 that 
would be expected for a test particle orbiting in the gravita- 
tional potential of the binary 



03 = - 



1 



where 



2r 2 i7kep dr \ dr 



(11) 



(12) 



is the phase-averaged non-Keplerian contribution to the grav- 
itational potential. Inserting the gravitational potential of the 
binary (equation|2]l in equations ( TT2l and dTTb yields 



16 \a/ \a 



n 



kep 



;4x 10 



-7/2 



fth 



(13) 



Figure [3] shows the final, quasi-steady-state disk eccentricity 
and ellipticity as a function of radius. For 1.5a < r < 5a, the 
quantity {e\ + e^} 1 ' 2 is roughly approximated by the function 

exp(-r/a). The eccentricity of the fluid in the central hole is executes quasi-periodic oscillations 



The measured, radially-averaged precession rate il prec agrees 
with the test-particle apsidal precession rate 03 at radius r = 
3.24 a, which falls within the range of radii where the disk 
appears to precess as a rigid body, and is close to the radius 
where the disk surface density peaks and most of the mass of 
the rigidly-precessing disk resides (cf. Figure [TJ. Therefore, 
the precession of the inner disk can be understood as arising 
from the non-Keplerian nature of the binary's gravitational 
potential. Exact agreement between £lp rec and 03 may not be 
expected at any specific radius because of the role of radial 
pressure gradients, which are neglected in the test-particle ap- 
proximation, that stabilize the rigid eccentric distortion. 

In addition to the linear trend ~ f! pre cf, the measured line of apsides 
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FIG. 4. — Time-averaged torque density dT jdr and the integrated torque 
T exerted by the binary on the disk at ~ 4000 binary revolutions. Time aver- 
aging was carried out over 45 binary revolutions, and reflects a quasi-steady 
state. The total torque experienced by the disk, 7"(oo), is positive. 

3.3. Torque Density 

The binary torques the disk by coupling gravitationally to 
nonaxisymmetric structures in the disk. In Figure @] we plot 
the time-averaged torque density 



dr \ a<p 



and the integrated torque 



, dT 



(14) 



(15) 



where, as before, <i> denotes the gravitational potential of the 
binary. The torque density can be of either sign. The strongest 
positive torque density is exerted at r < 2a. In this region the 
fluid density is relatively low, but the tidal force accelerates 
the fluid radially to highly supersonic velocities v r ~ vg ^> c s . 
A weaker, negative torque density is exerted at 2a < r < 2.4a; 
at larger radii the torque density oscillates around the zero 
point. 

To identify the origin of the oscillations in the torque den- 
sity, we work in the rotating frame in which the binary is at 
rest and carry out time averaging of the surface density in this 
frame. Figure [5] shows the surface density fluctuation 



ST,(r,e) = 



1 



(16) 



which selects the relative nonaxisymmetric component of the 
surface density distribution. The fluctuation is a two-armed 
tightly wound spiral density wave. Since the wave is station- 
ary in the rotating frame, its pattern angular speed in the iner- 
tial frame equals f\in- At large radii r 3> 2a, the wave obeys 
the Lin-Shu dispersion relation 



n 2 -m 2 (n- 



-n bin ) 2 +k 2 c 2 s : 



o 



(17) 



with ra = 2, where k is the radial wavenumber. 

The wave appears to be excited in the region a<r<2a and 
seems to refract at r ~ 2.5a where the surface density profile 
exhibits a jump. The oscillating torque density in Figure |4] 
directly follows from the stationary spiral pattern in Figure 
[5] Significant net torque is exerted only in the central, low- 
density hole and near the disk edge (r < 2a) where the spiral 
wave is nearly radial; the net contribution from larger radii is 
very small due to cancellation. 
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FIG. 5. — Time-averaged surface density fluctuation <5£ = (£ — (£))/{£) 
in the rotating frame in the which the binary is at rest. We have multiplied 
5S by (r/a) 1 for clarity. The binary point masses are located at r = ja and 
<f> = (0, n). The brightest regions are density maxima. The plot does not show 
the disk at radii r mm < r < 1.1a. 

The total, time-averaged torque experienced by the disk is 
positive, 

T(oo) «6x 1CT 4 GMaEo. (18) 

The binary-disk torque can be compared to the viscous torque 
communicated in the disk 



(19) 



-37ra I - I GMrY,. 



The two torques are of comparable magnitude when 
S/Eo ~ 0.64 a/r, which is roughly satisfied in the inner disk. 
This is of course expected if the disk has settled in a quasi- 
steady state characterized by a balance of viscous and gravi- 
tational torques. But the balance is not exact, because angular 
momentum is being transported by the (weak) radial accretion 
and density waves in the disk. 

3.4. Fluid inside the Central Hole 

Figure [6] shows that the fluid flow inside the central hole 
(r < 2a) is non-circular. Of course, in this figure we only see 
the high-pattern-speed m = 2 component, whereas the low- 
pattern-speed m = 1 component is stronger but not present af- 
ter time averaging in the rotating frame. Examining the in- 
stantaneous flow within the central hole, we observe pinball- 
like behavior; a fluid element located inside the hole receives 
a gravitational kick near the pericenter and is ejected outward 
on an eccentric orbit toward the inner edge of the disk. It even- 
tually collides with the inner edge at r ~ 2a. The momentum 
of the ejected fluid is transferred to the disk and a shock wave 
is driven into the disk. 

In § H] we interpret the torque transfer by periodic kicking 
of the eccentric fluid inside the central hole as a mean-motion 
resonance. In addition to the tightly-wound density wave ob- 
served in the frame rotating with the binary, a strong, one- 
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FIG. 6. — Same as Figure|5] with the addition of vectors representing the 
time-averaged velocity of the fluid in the rotating frame in which the binary 
is at rest. At r > a, the angular velocity approaches the negative angular 
velocity of the binary, dcfr/dt — > — f^bin ■ The brightest regions are density 
maxima. The plot does not show the disk at radii r m j„ < r < 1.1a, 

armed, long-wavelength density wave is driven at the apocen- 
ter of the inner edge of the disk. The latter wave precesses 
with a low pattern speed and hence disappears when time- 
veraging in the rotating frame. The long-wavelength wave 
appears to be associated with the propagation of the eccentric 
distortion in the disk. 

3.5. Accretion into the Central Hole 

The average mass accretion rate across the inner edge of the 
integration domain is 



MOmin) « 2.5 x 10- 4 (GMfl) I/2 So. 



(20) 



This can be compared to the nominal accretion rate in a disk 
without binary-induced torques 



d(r 2 n) 



dr 



dT v \ 



dr 



'6Tra[-) (GMr) I/2 £, 



(21) 



where we have assumed that in the disk without binary torque 
the surface density is a weak function of radius. Then we have 



M 

Mfree 



• 0.13 



-1/2 



(22) 



The accretion rate ratio in equation ( f22b should be evaluated 
at the outer edge of the zone where the binary forcing is signif- 
icant. This is the transition radius where the disk surface den- 
sity derivative changes sign, d(Y) /dr ~ 0; outside this radius, 
angular momentum transport in the disk is governed by non- 
gravitational, viscous torques only. Figures Q] and [4] suggest 
that the transition radius is at r ~ 3 a where (S) ~ 0.425 Eo- 
With this, M/Mftee ~ 0.18 < 1, which indicates that the ac- 
cretion rate into the central hole in the circumbinary disk is 
smaller than the accretion rate in a disk without binary torque 
with comparable density. 
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FIG. 7. — Time dependence of the accretion rate dM/dt = M — r J Vy'Edcp 
across the inner edge of the integration domain, v — r m j n = a, at ~ 4,000 
binary revolutions after the beginning of the simulation. 

The surface density profile of the disk exterior to the den- 
sity maximum at r ~ 3a is slightly steeper than the profile 
£ oc r~ l l 2 expected in a constant-Mf ree disk. However the sur- 
face density profile is significantly shallower than expected in 
a non-accreting, constant- r v j sc disk. In the absence of accre- 
tion t hrough the ce ntral hole, the disk behaves like a decretion 
disk (lPringlell991h . with surface density profile £ oc r~ l . This 
indicates that the outer disk has not yet settled in the viscous 
quasi-steady state; the simulated disk is neither exactly an ac- 
cretion nor a decretion disk. 

Figure [7] shows the time dependence of the accretion rate. 
This is the rate with which fluid flows across the inner bound- 
ary of the simulated domain; this flow should ultimately ac- 
crete onto either of the point masses or be otherwise removed 
(unbound) from the system. The accretion rate is quasi- 
periodic and exhibits notable low-frequency variation and pe- 
riodic outbursts. The peak accretion rate is a factor of ~ 4 
larger than the mean value. We discuss the implications of the 
time-dependent accretion rate in § [5] 

Figure [8] shows the Lomb-Scargle periodogram of the ac- 
cretion rate. The period ogra m was calculated u s ing th e defini- 
tions in [Scargli (QUI) and EornE&liiEnii CSl) 11 - The 
peak at ~ 2 x Slbin corresponds to the forcing frequency of 
the quadrupole of the tidal potential. The accretion rate time 
series is quasi-periodic as it can be decomposed into 10 com- 
ponents with evenly spaced frequencies uo ~ ^KQ\,i n , where 
K = 1,2, • ■ • , 10. The peaks at K = 3,4, and 5 are not signifi- 
cant. The most interesting is the fundamental frequency com- 
ponent at uj ~ ijilbin; it is responsible for the slow periodic 
modulation of the signal seen in Figure [7] 

3.6. Evolution of the Binary Separation 

The fluid accreting across the central hole transports angu- 
lar momentum at the rate 

L acc = 27rr 2 (v r V0l]). (23) 

We measure the average Lfl u ;d in the simulation and find 

iacc(r min ) » 3.3 x lO^GMaEo, (24) 

which is a half of the binary-disk torque (cf. equation [T8| ). 
Therefore the cross-hole accretion of (positive) angular mo- 
mentum is not sufficient to offset the loss of angular momen- 
tum to the disk, and the binary loses angular momentum on 
average. 

1 ' http://astro.uni-tuebingen.de/software/idl/aitlib/timing/scargle.html 
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FIG. 8. — The Lomb-Scargle periodogram of the rate of accretion across 
the inner boundary at r = r m ; n = a. The periodicity analysis was carried out 
over a 3 times longer time interval than shown in Figure[7j 

If the binary remains circular, evolution of the binary semi- 
major axis is governed by 



, _ (GMV^/a M 
Lacc l ~ 8 U M 



(25) 



where the minus sign in front of the torque on the left-hand 
side reflects T being defined as the torque exerted by the bi- 
nary on the disk. The rate at which the binary separation de- 
cays can then be estimated to be 



a 
a 



-0.003 



-7.3 x 10" 




(26) 



where M e( jge = 7r(2a) 2 I]o and M e d ge <C Mbh is assumed. The 
inferred binary orbital decay time is comparable to the vis- 
cous time at the inner edge of the disk, multiplied by the mass 
of the disk at the inner edge and divided by the binary mass. 
If the binary does not remain circular, the estimate based on 
our simulation is not strictly accurate, but the overall scal- 
ing should still apply. If the binary point masses are unequal, 
jMbh in equation (l26i l should be replaced by the reduced mass 
M X M 2 /M hia . 

4. DISCUSSION 

4. 1 . Torque Density from Tidal Forcing at m = 2 OLR 

The torque density profile in Figure|4]resembles the theoret- 
ical torque density profile for a density wave linearly excited 
at an OLR by a binar y's tidal potential in a circular, quasi- 
uniform-density disk dMever-Vernet & Sicardvl lT987). 12 To 
explore this similarity, we proceed to compare torques in the 
simulated disk with predictions of the linear response theory. 
In a circular, inviscid, non-self-gravitating, quasi-uniform- 
density Keplerian fluid disk that is forced by a rotating poten- 
tial multipole <£>„, of order m, torque density in the vicinity of 
the correspondi ng OLR can be written in term s of the Airy 
function Ai(jt) ( M eyer- Vernet & Sicardvll 19871) 



dT ^ T m / T m 
dr A™ V A„ 



(27) 



12 We are grateful to Steve Lubow for providing this reference, and for 
suggesting an outline for the present subsection. 



where r m is the radius of the resonance that is the solution of 
the equation ft(r m ) = [m/(m+ l)]Obi n and T„, is the approx- 
imate total torque exerted at an OLR in a pressureless disk 
dGoldreich & Tremaindl 1979b 



T m = mil 





dD 




s 




■'( 




dr 





d<S>„ 



dr Q-fl 



bin 



(28) 



Here, in a Keplerian potential, D = il 2 -m 2 (Sl- £7bi n ) 2 , and 
A m can be interpreted as the inverse radial wavenumber asso- 
ciated with the pressure wave excited at the resonance, 



3mr 2 QQ 



bin 



1/3' 



(29) 



The equal-mass binary analyzed here is special in that the 
potential component having m = 1 that corotates with the bi- 
nary is identically zero. Therefore, the outermost m = 1 OLR 
at r\ « 

2 2/3 

a is absent. The m = 2 potential component asso- 
ciated with the outermost OLR at r2 = (3/2) 2//3 a with orbital 
resonance flb m : fl = 3 : 2, in the special case of an equal-mass 
binary, equals 



<]), 



3_ 
16 



a 2 n 2 +o 



(:)" 



(30) 



Note that although r-i is the location of the resonance, the 
range of radii over which resonant forcing occurs is actually 
fairly broad, as is evident in Figure |H and is comparable to 
the radius itself. With this we obtain 



dT 49 a . a\. 

Ir- W 288" ^ r2 ^Y 2 M 



>~2- 



(31) 



where A 2 = 2- 2 / 3 (/!/r) 2 / 3 a. 

The approximate torque density derived from linear theory 
in equation OTT ) is only a crude approximation, because the 
simulated disk violates multiple premises on the which the 
derivation is based: 

1. The strength of the resonant torque associated with the 
m = 2 OLR is nonlinear, T2/[a 4 fl 2 tin Y,(r2)] > 1. In particular, 
nonlinear coupling with the eccentric m = 1 mode may modify 
the torque density. At the radius of the resonance the disk 
itself is intermittent, rather than steady. 

2. The derivation assumes that the forcing potential varies 
very slowly with radius, d In $2 /dr I/A2, but in the simu- 
lation, the forcing potential and the resonantly-excited wave 
vary on comparable length scales. 

3. The derivation assumes a quasi-uniform-density disk, 
dln(Y.)/dr <C I/A2, but in the simulation, the azimuthally- 
averaged density and the pressure vary more rapidly than the 
resonantly-excited wave. Since the resonance has a finite 
width, the resonant torque will depend on the density over 
a range of radii. 

4. Similarly, the linear derivation assumes an isothermal 
disk, whereas our disk is locally-isothermal and the sound 
speed varies as c s oc r" 1 / 2 . 

5. The derivation assumes a Keplerian potential, whereas 
the potential of the binary is stronger than Keplerian at small 
radii (e.g., see equation J4)). 

Each of these aspects can be, with some effort, incorpo- 
rated in an analytic derivation of the torque density. Instead, 
here we simply attempt to evaluate the torque density estimate 
according to equation OTT l after substituting the disk surface 
density (equation Q) evaluated at the resonant radius r 2 . The 



8 



CIRCUMBINARY MASSIVE BLACK HOLE ACCRETION DISK 



surface density equals S(ra) «6x 10 -4 So, and with this, the 
estimated torque density reaches maximum 



(32) 



at radius r « 1 .45 a. The peak torque density in equation d32i > 
is a factor of ~ 2.5 higher than the peak torque density in 
the simulation, which occurs at radius r w 1 .75 a (cf. Figure 
|4|. At larger radii, the theoretical torque density oscillates 
more rapidly than the simulated, that is, the wave excited at 
the resonance has a shorter wavelength (e.g., the fifth positive 
torque peak occurs at r » 2.87 a in linear theory, and at r w 
3.55 a in the simulation). This is not surprising, given that 
linear theory can severely underestimate the radial wavelength 
of the den sity wave excited at a n OLR when the forcing is 
nonlinear (I Yuan & Cassenl[l994l) . We can conclude that the 
torque density measured in the simulation is consistent with 
having arisen from a forcing by the binary's tidal potential 
at the m = 2 OLR. This implies that the same resonance can 
be responsible for disk truncation and maintenance of torque 
balance near the inner edge of the disk. 

4.2. Qualitative Interpretation of Disk Eccentricity 
Excitation 

We here provide a qualitative interpretation of the excita- 
tion and maintenance of disk eccentricity. Consider a disk 
fluid element orbiting the binary on a quasi-elliptical orbit 
with semimajor axis a and eccentricity e. The fluid element 
makes a pericenter passage at radius r p ~ (1 - e)a. Its angular 
frequency at the pericenter is approximately given by 



~ ^kep(fl) 



(1+e) 1 / 2 
(l-e) 3 /2- 



(33) 



16(l-e) 2 a 2 

When this frequency is close to the orbital frequency of the 
binary Slbin, the fluid element experiences a gravitational pull 
as it passes the pericenter point. Assuming that the eccen- 
tricity at the pericenter is e ~ 0.2 (see Fig. [3]), the resonance 
condition 

fi p (a, e l = ft bin (34) 

can be solved to obtain the numerical solution a « 1 .45 a. 

The fluid that arrives with a slightly larger a and slightly 
smaller e progresses in azimuth at a slightly lower rate than 
the binary point mass during pericenter passage. In this case, 
the binary imparts a gravitational "kick" to the fluid element. 
Transfer of momentum from the binary to the disk is evident 
in the positive torque density inside the central hole occupied 
by orbits close to the resonant condition in equation (l34l . as 
seen in Figure |4] The density wave (Figure ID is clearly ex- 
cited inside the hole. 

After departing the pericenter and proceeding to larger 
radii, the orbital frequency of the fluid decreases. In a time 
set by the relation 

/ (35) 

Jo 1 
the kicked fluid lags by a quarter of a phase behind the binary, 
and the sign of the torque changes, because now the other 
point mass is located closer to and just behind the kicked fluid. 
Therefore, a negative torque density is expected at larger radii, 
as seen in Figure |4] 

The point mass approximation employed thus far breaks at 
the apocenter. The simulation shows that the kicked fluid im- 
pacts the inner edge of the disk near the apocenter superson- 
ically which drives a shock wave into the disk. The kicked 



fluid imparts momentum to the disk which is then transported 
into the disk interior by a density wave. Therefore, the low- 
density fluid located in the central hole that experiences peri- 
odic kicks mediates the binary-disk torque coupling. 

After the impact at the apocenter, the kicked fluid loses or- 
bital phase coherence, and plunges back toward the pericenter 
as an eccentric stream. Consecutive kicks do not add up co- 
herently; the fluid does not retain the memory of its previous 
pericenter passages, which distinguishes it from t he collision- 
less ec centric mean motion resonance studied by iPan & Saril 
(2004). This puts into question the assumption at the heart of 
the standard treatment of binary-disk resonances that distur- 
bances are linear and that consecutive disk-planet passages 
can interfere coherently. Dispersal of the wave excited by 
each kick implies that the binary must resonate with the an- 
gular frequency at the pericenter alone, but not with any of 
the lower frequency components associated with the eccen- 
tric disk. 

The application of a positive torque near the pericenter and 
a negative torque near the apocenter acts to fortify the eccen- 
tricity of the inner disk. We speculate that this mechanism, in 
a linear form, is responsible for the transition from an ellipti- 
cal to an eccentric disk. Indeed, in an equal-mass binary, disk 
ellipticity can be linearly excited at the 2: 1 "eccentric" OLR 
(r ~ 1 .59a) by the mechanism de rived in iLubowl d!991l) and 
iD'Angelo. Lubow. & Bate] (120061) . 

Finally, we mention an alternative historical approach 
to disk truncation in circumbinary disks which is based 
on the behavior of the restricted three-body problem. 
The growth of orbital eccentricity of a collisionless third 
body orbiting a similar-mas s binary occurs in the vicin- 
ity of parametric resonanc es dHenonll 19651: ISzebehelvlll967l : 
Ruda k & PaczvnskilfT98ll) . If the orbit of the third body is 
interpreted as representing a streamline in the fluid circumbi- 
nary disk, the catastrophic self-intersection of resonantly- 
excited, e ccentric orbits leads to fluid removal and disk 
truncation dRudak & Paczvns1dll98lHBrwin & Sparkel[!999T: 
Pichardo, Sparke, & Aguilar 2005). We are not aware of an 
explicit demonstration that an equivalent of the three-body 
parametric resonances can truncate a fluid circumbinary disk. 

5. CIRCUMBINARY DISKS IN MASSIVE BLACK 
HOLES 

A massive binary black hole forms on a dynamical 
timescale at the conclusion of a galaxy merger, but can then 
persist at parsec or subparsec separations ove r a cosmologi- 
cal time (Begelman, Blandford. & Reeslll980l) . Initially, the 
binary decays by transferring orbital energy and angular mo- 
mentum to the stellar environment. Ho wever when t he separa- 
tion between the black holes becomes dPeterslll964t a circular 
binary is assumed) 



flgr(0~2500r g 



10 6 yr 



1/4 



(. 



M 



5 X 10 5 M r . 



1/4 



V 1O 6 M 
M 2 

5 x 10 5 Mr., 



-3/4 



1/4 



(36) 



where t denotes a time interval and r g = GM/c 2 with M = 
M1+M2, the emission of gravitational radiation drives the 
black holes to coalesce in time t . Note that the dependence on 
t and the total black hole mass is weak. In the gravitational- 
radiation-dominated phase, the binary will be very close to 
circular. 
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Consider a binary composed of low-mass massive black 
holes (M < 10 7 M Q ). Equation (f36t then implies that the bi- 
nary spends time t > 10 yr at separations a < a„ t ~ (10 3 - 
10 4 )r g . At such separations, the binary is sufficiently com- 
pact for it to be possible that a hot and ionized (T > 10 4 K) 
circumbinary accretion disk can encompass it. This type of 
accretion disk forms when interstellar gas accretes near a pre- 
existing binary. T he size of the disk is hmited by self-gravity 
at the outer edge dKolvkhalov & Svunvaevl [19801 : iGoodmanl 
I2003L and references therein). 

Since the circumbinary disk contains a central hole and 
the accretion across the hole and onto the individual black 
holes is reduced, the disk luminosity must be calculated self- 
consistently in any invocation of the Eddington limit. The 
standard "Eddington" accretion rate, calculated assuming a 
luminosity L ~ eMf lee c 2 , where e is a dimensionless efficiency 
factor, is of no relevance to the circumbinary disk. At binary 
separations a < 100r g , the binary decay timescale becomes 
shorter than the viscous r esponse time of the disk, and t he 
two components decouple (Milosavlievic & Phinnev 2005). 

Besides likely containing an unequal mass binary (Mi 
M2), the astrophysical hot circumbinary disk will differ in 
two significant aspects from the model disk in our simu- 
lation. The midplane temperature in the real disk at radii 
r ~ (100- 1000K will be in the range 10 4 - 10 6 K (e.g., 
iHubeny et al.ll200lh . implying a much cooler and thinner disk 
(h/r ~ 0.001-0.01) than in our simulation (h/r ~ 0.1). Fur- 
thermore, the locally isothermal equation of state assumed in 
our simulation will almost certainly be a poor approximation 
in the real disk; the inner edge of the disk will be heated by 
a combination of viscous torques and dissipation of waves 
driven by the binary. Nonetheless, the basic dynamical char- 
acteristics identified in the simulation should be present in the 
real disk; it is thus important to examine their implications 
for the detectability of compact massive black hole binaries 
embedded in circumbinary disks. 13 

The principal challenge to astronomical detection of an ac- 
tive galactic nucleus (AGN) containing a binary massive black 
hole is the confusion with normal AGNs that contain a sin- 
gle black hole. Our simulation suggests multiple forms of 
variability that can serve as indicators of a binary black hole 
in an AGN. The variability will be associated with: (1) the 
propagation and dissipation of density and shock waves in the 
bulk of the circumbinary disk; (2) supersonic gas dynamics 
in the central hole in the disk; (3) transfer of material across 
the hole onto the black holes; and (4) ultimate accretion of 
the material that has crossed the hole and has become cap- 
tured by one of the black holes onto the black hole, which 
may also proceed via a thin disk. Detailed modeling of the 
emission from these structures must be done to identify the 
most p romising indicators of b inary massive black hole AGNs 
(e.g.. lBogdanovic et al.ll2007l) . 

The amplitude of the variability produced by the strong su- 
personic shocks within the central hole can be estimated to be 
ihok ~ M ho i e v 2 r2bm, where M ho i e ~ a(h/r) 2 M edge is the mass 
of the gas inside the hole, v ~ (GM^/a) 1 ^ 2 is the orbital ve- 
locity, and as before M e d ee is the mass of the disk at the in- 
ner edge (we ignore dimensionless factors of order unity). 
This can be compared to the accretion luminosity of the disk 

13 An important caveat is our assumption that the disk is corotating with 
the binary. The dynamics of a counter-rotating disk will be very different and 
must be studied separately. 



idisk ~ GM hh M v i sc /a ~ a(h/ r) 2 (GM bh / a)M edge il. Evidently 
^hoie ~ ^disk, but the respective spectra will be considerably 
different; the emission inside the central hole will be harder 
th an that of the disk. 

lArtymowicz & Lubowl (1 19961) suggest that mass transfer 
from a circumbinary disk onto an eccentric binary will be 
punctuated by periodic outbursts. We, however, find that the 
mass transfer (see Fig. |7]l is quasi-periodic and highly vari- 
able even if the binary is circular; this is the case because 
the disk itself becomes eccentric. Temporal modulation of 
the cross-gap accretion in an eccen tric protoplanetary disk is 
simult aneously being reported by iD'Angelo. Lubow. & Bate! 
(2006), who simulate binaries with very unequal-mass com- 
ponents, Mo/Mi ~ 0.001-0.003. In a similar-mass binary, 
the principal period is half of the binary orbital period 

t ( t \ 3 ^ 8 / M \ ~ I//8 

f~ 23days U^r) [m^j 

x ( Ml \ ^ f Ml \ ^ (37) 
X V5xl0 5 M Q y V5xl0 5 M Q/ / ' 

but longer periodicities (e.g., ~ 4fbin) marked by outbursts of 
accretion across the hole can also be expected. The variability 
will be present on optical to X-ray wavelengths and can be 
detected by all-sky monitoring surveys such as LSST 

Accretion modulated by an e ccentric binary compan- 
ion (e.g.. ISillanpaa et alj [19881: iLehto & Valtonenl [19961: 
IValtaoia et alj|2000l) has been invoked as an explanation of 
the 12 yr optic al periodicity in the BL-Lacertae active galaxy 
OJ 287 (e.g. IS^paaetaflll996alrl iPursimo et alJ l200ol: 
IValtonen et al.l 12006). The optical light curve of OJ 287 is 
punctuated by tandem flares separated by ~ 1 1 yr of relative 
quiescence. The binary black hole models of OJ 287 involve 
a single or double crossing of one black hole's accretion disk 
by the other black hole; the semimajor axis of the black hole 
binary is assumed to be larger than the accretion disk. 

The quasi-periodic accretion across the central hole in our 
simulation (Fig. [7]) tentatively suggests an alternative explana- 
tion for the structure of the light curve of OJ 287. Evidently, 
a modulation of the accretion rate resulting in twin outbursts 
can take place in a circular binary surrounded by a coplanar 
disk. However a model of OJ 287 based on a circumbinary 
disk must explain the puzzling abse nce of a radio count erpart 
in the first of the two optical flares ( Valtaoja et al. 20Q3)- 

Finally we remark that if the disk is sufficiently massive 
while of course remaining gravitationally stable and geo- 
metrically thin, the torque imparted by the disk on the bi- 
nary can compete with stellar dynamical processes at extract- 
ing energy and angular mom entum from the binary and ac- 
celerating coalescence (e . g., [i vanov, Papaloizou, & Polnarev 
ll999tlGould & RixH "2000: Me rritt & Milosavljevicll2005L and 
references therein). At present we are unable to esti- 
mate t he expected mass of s uch a disk from first princi- 
ples. Armitage & Nataraian (2005) have shown that disk 
torques can excite eccentricity in a binary that has entered 
the gravitational-radiation driven inspiral. We speculate that 
the disk eccentricity that develops on a viscous timescale 
in our simulation ma y contribute to the excitation of the 
binary's eccentricity; Papaloizou, Nels on. & Massetl (1200 ll) 
have s hown that the two eccentricities are co upled. Re- 
cently, Bogdan ovic. Reyn olds. & Miller) (l2007h pointed out 
that cross-gap accretion will torque the angular momenta of 
the two black holes to align them with that of the disk, which 
may restrict the details of the ultimate gravitational-radiation 
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mediated coalescence to the regime not propitious to strong 
gravitational recoil. This suggests that the consequences 
of geodetic precession that is present if the black holes are 
rapidly rotating should be explored in future extensions of our 
study. 

6. CONCLUSIONS 

We carried out a grid-based, two-dimensional hydrody- 
namic simulation of a thin, locally-isothermal, corotating, 
non-self-gravitating, viscous accretion disk around a Newto- 
nian equal-mass binary on a fixed circular orbit. The disk is 
evolved over multiple viscous times of the inner disk. We 
study the quasi-steady-state structure of the circumbinary 
fluid flow and the mechanics of the binary-disk torque cou- 
pling. 

We here summarize the main conclusions of this work. 

1. Fluid density within twice the binary's semimajor axis 
is significantly reduced. The point masses orbit within the 
central low-density hole in the disk. The disk first becomes 
elliptical and then eccentric. The eccentricity of the disk is 
strong and precesses slowly. Radial fluid motion within the 
central hole is supersonic and time-dependent. Strong shocks 
are driven into the low-density fluid by the tidal field. 

2. The binary exerts torque on the disk by imparting gravi- 
tational kicks to the eccentric low-density fluid inside the cen- 
tral hole. Angular momentum is transferred to the disk as the 
kicked fluid impacts the edge of the hole. A two-armed spiral 
density wave that is stationary in the rotating frame in which 
the binary is at rest is driven into the disk. The structure of 
the two-armed density wave and that of the binary torque dis- 
tribution in the disk suggest that the binary forces the disk at 
the m = 2 outer Lindblad resonance, where the orbital period 
of the disk is 3/2 times the binary's orbital period. 

3. The accretion into the hole toward the two point masses 
is reduced relative to the rate expected for a disk with the same 
surface density and a torque-free inner edge. However, sig- 
nificant mass transfer from the disk onto the black holes does 
take place. The mass transfer rate is quasi-periodic and punc- 
tuated by outbursts. The accreting fluid carries angular mo- 



mentum, but the amount of angular momentum accreted onto 
the binary is insufficient to offset the loss of angular momen- 
tum to the gravitational torque. If the binary remains circular, 
its separation decays on a viscous time of the disk reduced by 
the ratio of the mass at the inner edge of the disk to the binary 
mass. 

4. The quasi-periodic accretion across the hole, the shocks 
in the low-density fluid inside, and the disk eccentricity are 
sources of variability that can be used to distinguish an AGN 
containing a binary massive black hole from a regular AGN 
containing a single black hole. Detailed modeling of the emis- 
sion from these structures is required to identify the most 
promising indicators of binary massive black hole AGNs. 
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